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I. INTRODUCTION 

Recently, there has been a revival of the vortex picture of confinement Q]. Lattice calcu- 
ations performed in the so-called maximal center gauge with subsequent center projection 
^ (which provides a simple method for identifying the center vortices in Yang-Mills con- 
figurations) have given substantial support to the idea that center vortices are the relevant 
infrared degrees of freedom responsible for confinement. The deconfinement phase transi- 
tion in the vortex picture takes the guise of a percolation transition 0, |^. Furthermore, 
the topological charge of the vortices can be understood in terms of their self-intersection 



number 



2| and the spontaneous breaking of chiral symmetry can be generated 
through vortex effects III | ■ 
On the other hand, while a successful phenomenological model of vortex dynamics has been 
developed investigations of the dynamical properties of Yan g-M ills vortices have 



only recently been launched, cf. 



LS oi tne aynamicai properties oi lan g-iy i 
1^ 13 1 for thick center vortices and jlJ] 



for thin vortices 



of arbitrary flux. In the present paper, we expand on these investigations, which only 
considered gluonic fluctuations to one-loop order, by evaluating the one-loop energy density 
of an infinitely thin static magnetic vortex including quark fluctuations. For this purpose 
we use the Schrodinger picture. We solve the Schrodinger equation for gluonic fluctuations 
and fermions moving in the vortex background field. From the obtained eigenenergies we 
calculate the static energy density of the vortex as a function of its magnetic flux. We will 
show that the inclusion of fermions substantially changes the energy density of the vortex. 
In particular we will find that the fermions raise the center vortex effective energy compared 
with the perturbative vacuum. 

The organization of the paper is as follows: In the next section we study the Schrodinger 
equation for the gluonic fluctuations around the vortex background field. In subsection III Al 
we study the Hamiltonian of the gluonic fluctuations in the vortex background and discuss its 
self-adjoint extensions in subsection III Bl Subsection III (^l is devoted to the incorporation of 
Gauss' law and to the selection of the appropriate self-adjoint extension of the Hamiltonian 
of the gluonic fluctuations. In subsection III Dl the gluonic contribution to the energy density 
of the vortex is presented. 

In section irm we study fermions in the vortex background field and calculate their contribu- 
tion to the energy density. The physical interpretation of our results and some concluding 
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remarks are given in section HVl 



II. GLUONIC FLUCTUATIONS AROUND THE VORTEX BACKGROUND 
FIELD 

A. The Hamiltonian of the gluonic fluctuations 

Our aim is to evaluate the one-loop corrections to the energy density of a static SU(2) 
vortex background field. We consider an infinitely thin vortex of magnetic fiux $ on the 
3-axis described by a gauge potential 

Mp^v) = ^—e^ . (1) 

^ UP 

where the polar coordinates (p, (f ) parameterize the (1, 2)-plane in space. The normalization 
has been chosen such that odd integer values of $ describe a center vortex, i.e. any Wilson 
loop W in the (1, 2)-plane encircling the origin takes the value W = —1. The magnetic field 
associated with ((T)) is 

B\x) = 8''^'^-^5{x,)8{x,)e, . (2) 

Our starting point for the evaluation of the one-loop energy density is the Yang-Mills Hamil- 
tonian in Weyl gauge {Aq = 0) in the presence of an external current 



H= I d^x 



^E^E^ + ^B^Bt-j^A} . (3) 
By decomposing the gauge field Af into the vortex background field Af and fiuctuations rjf, 

At = At + vt , (4) 

dropping terms of higher order than quadratic in the fiuctuations rj (one-loop approxima- 
tion), adjusting the external current such as to induce the background field A"" at the classical 
level, i.e. 

DfF^ = -j^, Df = 6''% + gr"'A'r (5) 
and finally choosing the background gauge 

Dfv^ = (6) 



for the fluctuations, one obtains tlie Hamiltonian (in the "coordinate" representation) 



H 



i/cl ^ ^(1) ^j^J^ 



X 



X 



BfB^ - 2i^Af\ and 



(7) 



52 



_|_ ^a/\ab^b 
' li ij Ij 



Here the operator O reads 



k 5 



(8) 



where {D'^Y^ = D'^'^Df. In this approximation, the ground state wave functional takes the 
Gaussian form 

^ri] = exp (-^ y VKV^ (9) 
with the spectral representation for the covariance K 



N 

Here (p{N) denotes the eigenfunctions of O, cf. (jH}, 

d(j){N) = 0{N)(j){N), 



(10) 



(11) 



where denotes a (not necessarily enumerable) complete set of quantum numbers. Fur- 
thermore, the eigenfunctions 0(A^) have to fulfill the "background gauge" condition 



D(f){N) = 



such that the wave functional fulfills Gauss' law, 



(12) 



(13) 



to the desired order, namely 0{g^), and is an eigenfunctional of the Hamiltonian ((7j) with 
eigenvalue 



E = E"' + ^ / d^x 



X] 



N 



d X e(x) 



(14) 
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where e(x) denotes the energy density. To calculate e(x), the set of eigenfunctions must 
be specified. As a first step, it is useful to reduce the problem by decomposing the functions 
(j) (and the quantum numbers N) as 

0^(A^, x) = Uj{i)e^{a)fMixi,X2) exp(iA;x3) exp(imy9) , A; G R, m G Z (15) 
with polar color and Lorentz unit vectors 

e(l) = (l/v^,i/v^,0), e(2) = (l/v^, -1/^2,0), e(3) = (0, 0, 1) , 

(16) 

n(l) = (1/v^, i/v^, 0) , u{2) = (1/v^, -i/v^, 0) , u{3) = (0,0,1) , 
in terms of which the eigenvalue equation (|TT|) reduces to 

I3{a,i,m)fM = B{M,a,i,m)fM (17) 
where the operator B is defined by 

A/ . N 1 d d (m — A(a)$)^ x„ r/ x <-/ x /,^\ 

B(a, I, m) = — —p— + !: - 2X(t)X(a)2TT^6(xi)6(x2) (18) 

p dp dp p^ 

and its eigenvalues B{M, a, i, m) are related to the ones of O, cf. (jHI), by 

0{N) = B{M,a,i,m) + k\ (19) 

Note that the color quantum numbers a and the Lorentz quantum numbers i enter via the 
eigenvalues A corresponding to the vectors (fTI)|) . namely 

A(l) = 1 A(2) = -1 A(3) = . (20) 

The additional constraint on the /m imposed by Gauss' law, eq. (fT3|l . will be dealt with 
later. As the operator B, cf. (I18|) . is not well defined on the z-axis, we restrict space to 
]R'^\{2; — axis} and take the delta functions into account by adopting appropriate boundary 
conditions. 

At this stage, it is useful to introduce a scale M which has the dimension of a mass. This 
allows us to make eigenvalue equations dimensionless and, when we determine boundary 
conditions, we make an expansion of functions for distances to the 2;-axis small compared to 
the scale 1/M. With this scale, we define a new variable x := pM and obtain a dimensionless 
operator 



For a complete definition of tlie set of functions /m, a set of boundary conditions in 
L^(R>o,a;dx) must be specified for a; — > 0, such that C is self-adjoint. In order to sys- 
tematically construct all self-adjoint extensions, we initially restrict the domain D{C) to 
D{C) = C^(]R>o) on which C* is a well defined symmetric operator, i.e. 

{f\C9) = {Cf\9) forall/,^eD((7) . (22) 

Self-adjointness, by contrast, means more; namely, that C is equal to its adjoint C"^ , defined 
as follows. Let /^(C*^) be the set of functions g for which h G L^(]R>o, xdx) exists such that 

{9\Cf) = (h\f) for all / G D{C) . (23) 

Then, for each g G D{C'^), one defines C"^ g = h. As D{C) = C^(]R>o), one sees that in 
the sense of weak derivatives ~^-^x-^g is locally in L'^iH^Q, xdx) and therefore also locally 
in L^(]R>o, xdx). From this one can conclude that {xg')' is locally in L^(]R>o,dx) with the 
usual measure and, with Sobolev's lemma that g' is absolute continuous. This allows 
an integration by parts in (j^Hj) which gives 



C^ = -^±x— + - (24) 
X dx dx x^ 

D{C'') = {/ G L^(]R>o, xdx) I/' absolute continuous, e L^(]R>o, xdx)} , (25) 

where the derivatives are now ordinary derivatives and we introduced 5 = \m — A(a)$|. 
Note that D{C) C D{C'^) as it should be for the symmetric operator C . In the following, 
we construct all self-adjoint extensions of C . In the cases where we obtain more than one 
self-adjoint extension, which corresponds to different possible boundary conditions on the 
z-axis, the choice of a definite self-adjoint extension is a question of the physical situation 
one wishes to describe; mathematically, all choices are equally valid. 



B. Self-adjoint extensions 

To determine self-adjoint extensions of C*, it is necessary to consider the deficiency spaces 7i± 
of C*, i.e. the spaces of zero modes of the operators C*^ =F i, with their respective dimensions 
(deficiency indices) n± [1^]. The operator C is called essentially self-adjoint if n+ = n„ = 0, 
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thus admitting a unique self-adjoint extension given by its closure. With / G D{C^), cf. (j2SI), 
and 



C^f = ±if (26) 



one can conclude that / G C^(]R>o). Two linearly independent solutions of C^f = if are 
Is (^"^(1 ~ (^'^i^ ~ '^)^) ' where Ks and Is denote the modified Bessel functions. 

With the condition / G D{C'' ) only the solution 

,72 



f4x) = Ks{^il~t)x] (27) 



for < 6 < 1 remains. Therefore 



n+ = {c/+|cG C} for < 5 < 1 
= dl for (5 > 1 . 



(28) 



Correspondingly, we obtain 



H- = {c/_|c G C} for < 5 < 1 
= dl for 5 > 1 



(29) 



f.{x) = Ks{^il + i)x]. (30) 



with 

f (rr\ = kJ 

As already discussed above, this means that there exists a unique self-adjoint extension for 
6>1. 

For < 6 < 1 the deficiency indices are equal: = n_ = 1, such that there is a one parame- 
ter family of self-adjoint extensions Ca {a G (— vr, vr]) which are in one-to-one correspondence 
to isometrics Ua '■ f+ ^ exp(iQ;)/_ from to in the following way: 

D{C^) = {/ + c(/+ + UJ+)\f G D(C7), c G C} 
a(/ + c(/+ + f/„/+)) = (7/ + c(i/+ - iUJ+) . 

In Appendix it is shown that the functions / G D{C) in the closure of C are absolute 
continuous close to the origin and fulfill lim^^^o f{x) = 0. As the function /+ is singular at 
the origin, the boundary conditions at x = is determined by the behavior of this function. 
In the following, we determine these boundary conditions for the cases < 6 < 1 and 6 = 
separately. 



< 5 = \m- A(a)| < 1 

Expanding / G D{Ca) for small x, one obtains 



f{x) = A{x^ + /3{a)x-^) + o{x^) , as a; ^ (31) 
for a constant A E <C and with the one-to-one map 

^^"^ • Ti-5) cos(f + f 5) • 

Since CJs{k,x) = h?Js{k,x) and CNs{k,x) = k,'^Ns{k,x), one can make the following 
ansatz for the scattering states, 

AJs{kx) + BNs{kx) 

and again expand for small x. Comparing this with (j31|) yields the condition 

A cos(f + f (5) 



B K25cos(f - f5)sin(57r) 



cot(57r) . 



As bound state, only the solution Ks is admissible, since CKs{kx) = —k'^Ks{kx), and 
Ks is square-integrable. If one expands Ks for small x and again compares with (|31|). 
one notices that Ks with 

'cos(f + f5)' 



cos(f - U) 



is a possible eigenf unction, if 



cos(f - f 5) 
is satisfied. 

5 = \m-\{a)\ = 

Expanding / G D{Ca) for small x, one obtains 



f{x) = A (in (^1) + + o(x°) , as X ^ (32) 

for a constant A G C and with the one-to-one map 

/5(a) := 7 - - tan 
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where 7 is Euler's constant. For the scattering states, one makes, as above, the ansatz 



AJo{kx) + BNq{kx) 

which, after expanding for small x and comparing with ()32|). yields the condition 

A 2, , , 1 /( 

— = m K tan - 

B TT ^ ' 2 v. 

The possible bound state is 



Kq{kx) 

provided that 



4 V2 

is satisfied; correspondingly, one has the eigenvalue 



K = exp ( — Y tan ( — ) ) (33) 



—K^ = — exp I tan , 

^ V 2 V2 

For a = 71 there is no bound state and the solutions are regular at the origin {B = 0). 



TT /a 



C. Gauss' law and choice of self-adjoint extension 

In the previous section, for the cases < 6 = \m — A(a)$| < 1 a whole family of self-adjoint 
extensions of the operator C was obtained, labeled by one continuous parameter a; each 
value corresponds to a different choice of boundary conditions reflecting different physical 
situations. In this section, a definite choice for the parameter a will be motivated. 
The operator B, cf. |T8|) . can be interpreted physically as the Hamiltonian of a par- 
ticle with effective angular momentum (grand spin) 6 = m — A(a)<l' in a potential 
—2X{i)X{a)2TT^6{xi)6{x2). Therefore, a bound state can only arise if the potential term 
is negative, i.e. X{i)X{a) > for positive flux $ > 0. This is the case for X{i) = X{a) = 1 
and for X{i) = X{a) = — 1. It will now be demonstrated that Gauss' law forbids bound states 
even in these two cases. 

In order to fulfill Gauss' law, it is necessary to restrict the set of functions 0(A^) in each 
subspace of definite eigenvalue 0{N) to those linear combinations which satisfy eq. (fT2|l . 
Since the corresponding operator Df^ does not mix functions 0(A^) with a different color 
quantum number a, such linear combinations can be constructed for each a separately. The 
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same is true for the momentum in 3-direction k and the angular momentum quantum number 
m. Concentrating on the bound state solutions in the case A(a) = X{i) = 1, one has, using 
the decomposition (P3|). 

Df(j)%N) = 0^{pd,~m + $)/m = (34) 
which for the bound state solutions /a/(p) = K\m~<s>\{i^p) would require 

{\m - ^ - m + ^)K\ra-<i,\{Kp) - KpK\m-<i>\+l{Kp) = 0. 

For this to be true, both coefficients must vanish; in particular, k = 0, which contradicts 
the condition k > 0. Therefore, there can be no bound state for A(a) = X{i) = 1. The 
same is true in the case A(a) = X{i) = —1. For $ G Z, this singles out a definite choice 
of the self-adjoint extension, namely a = tt, since this is the only choice which yields no 
bound state. In this case, B = 0, i.e. the functions are regular at the origin. It should be 
noted that for non-integer fiux $, in general there are several possibilities with no bound 
states. In the following, however, regularity at the origin {B = 0) will be demanded even 
for non-integer fiux. 

From the physical point of view, this outcome is natural. Gauss' law enforces gauge in- 
variance; gauge invariance in particular implies that infinitely thin integer fiuxes, which 
are pure gauges in the theory without quark degrees of freedom, should be unobservable. 
Therefore, the gluonic excitation spectrum in the presence of such a background field should 
be equivalent to the free case. However, it should be emphasized that other interpretations, 
corresponding to other choices of self-adjoint extension, are possible and should not a priori 
be dismissed as invalid; e.g. the thin vortex fiuxes treated here may be used as idealizations 
of physical, thick fiuxes which are not pure gauges. Then, a different choice of self-adjoint 
extension may be useful to mimic the effects of vortex thickness. 



D. Energy density of a vortex to order 

For the particular choice of self-adjoint extension adopted above, in the sectors with A(a) 
±1 the energy density associated with a vortex differs from the one in the vacuum by 



e± = 2— I I dk ^ y^(/tM)2 + P (j|„zp$|(Kx))^ - (J|^|(/tx))^ , (35) 



27r 
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(for X > 0) where the prefactor 2 comes from the two physical polarizations allowed by the 
Gauss'law constraint. In the sector with A(a) = 0, the energy density (compared with the 
one of the vacuum) obviously vanishes. 

Let us first consider an integer-valued flux $ G No, A(a) ± 1. For $ G No, one has 

oo oo 

^ (J|„l(Ka;))^ = 1 and ^ {J\^^^{Kx)f = I . (36) 

m=— oo m=— oo 

Thus, the energy density vanishes for all x > and therefore 

e(p) =0 for p > . (37) 

The invisibility of thin integer vortex fluxes $, which represent pure gauge backgrounds, 
thus becomes manifest. 

Let us now consider fluxes $ G (0, 1). Consider first the energy density (j35|) in the sector 
given by A (a) = 1: Rewriting the square root in eq. by a proper-time integral 



oo 

i/a = — lim / d 

Smin^O J 



as 



(38) 



m=—oo 



where Smin is an (ultraviolet) proper-time cut-off we find 

oo oo 

ds 27r 



{j\m-<S>\{l^x)y - (J|m|(ra)) 

The integrals over k and k can be performed by using the relation 

d« e-^(*")V,(«x) 



-s({A?K)2+fc2) 



X 



K 



1 ^2 f \ 

^ = ^e2sA?2 J ( ^ I 1/ = Ittj — $1 or 1/ = \m\ , (39) 



cf. eq. 6.633.2 in yielding 

oo 

e^ = -2M^ / ds(7rs)-i/2A^ 

ds 27r V s ^ 27r 2sM2 



1 1 



m=— oo 



|m— 



\2sM^ J 



/|ml I - 



X 



(40) 



By substituting t 



e-, 



= f^, one arrives at 



? /2tdl/2t^lt 

m=— oo 



27r2o4 



tmaa: / OO \ 

y" dtt^/^^t^^/^e"* |^/,i,(t) - /o(t) + [/„^-$(t) + Im+^{t) - 2/„(t)] j 
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oo 



After inserting the integral representation jioj l 



into e+, the first terms of the integral representations of the different modified Bessel func- 
tions cancel, so that the only non- vanishing contribution reads 

tmax / 



r2^4 



oo 



j dtti/2At3/2g-tl Ij c^2e-*™^^(")sin($7r)e- 



oo °? 

+ J2 d^e-*™'^(^) (-cos(m7r)sin(<l>7r)e-(™-*)^ + cos(m7r)sm(<l>7r)e-(™+*)^) 



m=l 







2 sin(<l>7r) /■ ^^^1/2^^3/2 /"^^g-t(cosh(.)+i)COsh(2;[2-$]) 



27r2p^ 7 rft 7 cosh(f) 



Partially integrating with respect to t and noting that the surface term vanishes in the limit 



i^max oo, one obtains 



2 sin(<l>7r) /■ /^^g-t{cosh(.)+i)COsh(2;[2-$]) 



^ (27rp2)2 ^ J J cosh(f) 



For p > the limit Smin — *■ implies tmax oo. Then the integration over t can be carried 
out with a finite result, 



oo 

2 sin(<l>7r) f ^ 1 cosh(z[| - $]) 



(27rp2)2 J (Cosh(z) + 1)2 COsh(f ) 





oo 



2 sin(<l>7r) , f , e^(V2-*) + g-^^i/^-*) 
^ -4 / dz 



(27rp2)2 7 (e^/2 + e-^/2)(e^ + e-^ + 2)2 ■ 



Substituting x = yields 

oo 

2 sin($7r) f dx x^/^-* + a;i>-i/2 



1 

00 



2 sin(<l>7r) f , x^-* + 
' da; 



(7rp2)' TT 7 (X+1)5 

1 

2 sin(<l>7r) r(3 - $)r(2 + $) 

(VF ^ r(5) 

^ ^ -$(1-<I>2)(2-<1>) for $g(0,1). 



p4 127r2 
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FIG. 1: Energy density eges as a function of magnetic flux (/>. 



For e_, i.e. A(a) = —1, one obtains the same result in a completely analogous fashion. 
Furthermore, according to ()35|) the energy density is invariant under integer shifts in $. 
The total energy density 



eges = e+ + e_ = -— 0(1- 02)(2-0) for $ = iV + 0, iV e No, G [0, 1) (41) 



is depicted in Fig. The fact that the energy density is proportional to is natural on 
dimensional grounds. Furthermore, for integer-valued the fluctuations are equivalent to 
the ones in the vacuum, as expected for background fields representing pure gauges in the 
theory without quarks, i.e. in the presence of fluctuating fields in the adjoint representation 
of the gauge group only. Such fluctuations are not sensitive to a thin integer flux vortex. 

III. FERMIONS IN THE VORTEX BACKGROUND FIELD 

In this section we calculate the energy density of the quark fluctuations in a vortex back- 
ground field. The Hamiltonian of the quarks tplx) in the gauge field v4^(x) satisfying the 
Weyl gauge Ao{x) = is given by 



This result agrees with the one obtained in ref. [1J|. However, there the freedom in the choice of the 
self-adjoint extension was not discussed and a particular choice was made ad hoc. 



1 1 




(42) 
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where 



h{x) = a{p + gA) + f3m (43) 

is the Dirac Hamihonian with ol,(3 denoting the usual Dirac matrices. In the following 
we will, for simplicity, consider massless quarks (m = 0), which is a good approximation for 
light quark flavours. Furthermore, we will use the Weyl representation of the Dirac matrices, 
such that the Dirac operator, eq. (jlBI), becomes block diagonal, 

fc= ( " ° . D = ^{p + gA) (44) 

where cr are the usual Pauli spin matrices. For the gauge potential of the vortex field defined 
by eq. the operator is diagonal in color space 



Dab = SabD{a), D{a) = tr (^-iV + ^^e^) (45) 

where A (a) = 1,-1 for a = 1,2. Due to the geometry of the vortex field, it is convenient to 
use cylindrical coordinates {p,{p,z) in which the operator, eq. (jlSI), reads 



-id, ie-'^ (-a.- i (/ + 

where / = —id^ is the orbital angular momentum. It is seen that the vortex field generates 
an intrinsic angular momentum (color spin) of the quarks 



s = h{a)<!>. (47) 

For center vortices, where $ = 1, and, since |A(a)| = 1, the total (grand) angular momentum 
of the fermions J = l + s becomes half- integer valued, which is reminiscent of a boson- fermion 
transmutation. 

To find the spinor eigenf unctions u of the operator D{a), cf. fl46|) . we make the following 

ansatz. 
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u{a,k,m,N';x) =e'''' \ ^ ' , A; G R, m G Z (48) 

e'(™+i)^?/'(a,A;,m,iV';p) 

where the quantum number N' will be specified further below. With this ansatz, the eigen- 
value equation for the operator -D(a), eq. pUj) . reduces to 



^ ' P M , , , \ = E{a,k,m,N') {^^ ' (49) 

-i(9,-p -k J \^{a,k,m,N'-p)J \^{a,k,m,N';p)^ 

where we have introduced the quantum number 

z/ := m + ^ ' (50) 

of the grand angular momentum J = Z + s. Note that the eigenvalue equation ()49p is 
equivalent to the Dirac equation in two space-time dimensions for a particle with mass k. 
In the standard fashion, eq. (jl^ can be reduced to a second order differential equation for 
the radial functions Xi which takes the form of Bessel's differential equation. The solutions 
of the eigenvalue eq. PH]) for E'^ > k"^ (scattering states) are given by 



1 X 1 


1 = AT ( 








1 = AT j 







VETk{erJ,,{Kp) 
iyE^(e)-+V,(,+i)(«:p) 



W-E - k{eTJ,,{np) 
^_^ + A;(e)"+iJ,(,+i)(/€p) 



for E > 



for ^ < 



(51) 



where A/" is a normalization constant and k = ^/E"^ — k^ is the wave number in the plane 
perpendicular to the center vortex. Furthermore, Ja{x) denotes the ordinary Bessel functions 
and we have introduced the parameter e = ±1. The wave functions, eq. ()5ip . are solutions 



of eq. ()49p for both values of this parameter and for the eigenvalues E = ±\^k?+W. 
For u = n E X the solutions with e = ±1 are linearly dependent, since the Bessel functions 
satisfy the relations J_„(x) = (— 1)" J„,(x); thus, one should take Bessel functions and 
Neumann functions as linearly independent solutions. For z/ ^ which holds in particular 
for a center vortex field, for which $ = 1, the solutions with e = ±1 are linearly independent, 
i.e. there are two linearly independent solutions for each eigenvalue. However, requiring the 
functions to be square integrable on the 2;-axis selects out a single solution (that is, either 
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e = lore = — 1) except for — 1 < z/ < 0. The latter interval includes the case ^ = — | which 
occurs for a center vortex, $ = 1 (when A (a) = 1 and m = — 1 or when A (a) = —1 and 
m = 0). In this case, both solutions e = ±1 are square integrable on the z-axis. 
In the following, we will investigate the case — 1 < z/ < in more detail. The point is that 
we have to impose appropriate boundary conditions on the z-axis, where the vortex field 
has a singularity, so that a unique solution survives. The boundary conditions have to be 
chosen in such a way that the operator in eq. ()45|) becomes self-adjoint. The procedure 
works analogously to the gluonic case: One restricts the domain of definition of the operator 
D{a), eq. ()45|1 . in a suitable way, such that it becomes symmetric. Then one determines 
the deficiency spaces of this operator and, from this, all its self-adjoint extensions. The 



construction of all self-adjoint extensions of the operator in eq. ()49|1 was given in 
The extensions are characterized by a parameter 6 G [0, 27r) and defined by the following 
boundary conditions on the eigenfunctions. 



\^(A?p)-''-icos(J + f)y 

where M is a mass scale introduced in order to have dimensionless quantities like Mp. The 
requirement of self-adjointness restricts the choice of the boundary conditions on the vortex 
axis {z-a.xis). For the boundary conditions ()52|1 . a unique solution follows for each value of 
6 G [0, 27t), which in the case of > is given by 



= ^[1 + (-I)'" sin(2p) cos(z/7r)]-i/2 

v/ETI[sin(/i)J,(/cp) + (-!)- cos(p)J_,(ft:p)] \ ^^^^ 
i^/E^[sm{^i)J,+l{Kp) + (-!)"+! cos(p)J_(,+i)(/s:p)] 



where the quantity p is related to 6 by 



tan - + - = (-1)'" — — ^ ^tan(u). (54) 

V4 2 J ^ ' \E-k) \2k) r(z/ + l) ^ ' 

In addition to these scattering states, for | < 6* < ^ bound states, i.e. square integrable 
eigenstates exist jl?!. Iisj]. 



16 



Finally, let us remark that for u < —1 and for z/ > 0, where the parameter e = ±1 is already 
uniquely determined by requiring square integrability of the wave function on the z-axis, the 
application of the above procedure for finding self-adjoint extensions of the Dirac operator 
is trivial. In fact, one finds that the deficiency spaces of this operator are empty and that 
hence this operator is essentially self-adjoint. 

In those cases where there is more than one self-adjoint extension, one has to select one 
on physical grounds. In reference it was shown that only for ^ = ±| a vortex on the 
z-axis can be consistently described in the sense that the eigenf unctions of the self-adjoint 
extensions of the Dirac operator and its square fulfill the same boundary conditions. It is 
still open which sign of 6 has to be chosen for each flux $. The sign of 6 can be fixed by 
requiring that the energy of the quarks in a vortex field should be a periodic function in 
the flux $ of period 2. This is because fermions should not feel a Dirac string, $ = 2. It 
turns out that this requires choosing 9 such that the assignments listed in the table below, 
cf. ()56|) . result. Let us emphasize that this implies different signs for different ranges of flux. 
The same self-adjoint extension of the Dirac operator was obtained in Q by demanding 
that the expectation value of the resulting baryon number 



(55) 



properly transforms under charge conjugation, (A^) —{N), and is a periodic function in 
the flux $ with period 2. 

It turns out that this choice of the parameter 6 corresponds to choosing the minimal irreg- 
ularity of the eigenfunctions on the vortex axis. 





< $ < 1 


1 < $ < 2 


A(a) = 1 


e = 1 for m = 0, 1, ... 

e = —1 for m = —1, —2, ... 


e = 1 for m = —1, 0, ... 
e = —1 for m = —2, —3, ... 


A(a) = -1 


e = 1 for m = 0, 1, ... 

e = —1 for m = —1, —2, ... 


e = 1 for m = 1, 2, ... 
e = — 1 for 171 = 0, —1, ... 



(56) 



Furthermore, the same choice of self-adjoint extension was obtained in reference j2l| by 
choosing boundary conditions of the Atiyah-Patodi-Singer type for a finite radius of the 
vortex flux and subsequently letting this radius go to zero. 
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A. The energy density of quarks in a vortex background 

In the following, we calculate the energy density of fermion fluctuations in the vortex gauge 
field. The energy density is defined by 



e{x)=J2 Yl E{a,N,n){'^''y{a,N,n;x)'$''{a,N,n;x) (57) 

"=1.2 E{a,N,n)<0 

where E and \1' denote the energy eigenvalues and eigenf unctions of the Dirac Hamiltonian, 
eq. ()43p. Due to the block structure of the Dirac Hamiltonian, these eigenvalues are given by 
both the positive and negative eigenvalues E appearing in eq. ()49p. and the energy density 
becomes 



^ ( V E{a,N)u\a,N;x)u{a,N;x) 

a=l,2 ^ ,r Jr^.r^ 



:e{a,E<0\x) 



(58) 

- ^ E{a,N)u''{a,N;x)u{a,N]x) 

N,E{a,N)>0 



::e(a,£;>0;a;) 



with u denoting the wave functions corresponding to the energy eigenvalues E. We have 
assumed here that the wave functions are properly normalized (as scattering states), 

J pdp v) (a, K, x) u (a, A^, k', x) = — 5 — . (59) 



This fixes the normalization constant to 



We have to choose the eigenf unctions such that they respect the boundary conditions cor- 
responding to the minimal irregularity on the vortex axis (j56|) . 

With this choice of the eigenfunctions, one finds after a straightforward calculation the 
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energy density for fluxes < $ < 1, 



/ oo \ 

ei{^;x) = J (^k j KdK Vn^ + k^h ^ Jin+^/2\M + J-^/2M - ^1/2^] 

-oo ^ ™=-°° / 



while for fluxes 1 < $ < 2 one obtains 



(61) 



/ oo 

e2($;x) = -— dk ndK Vn^ + k^ h ^ Jfm+^/2\M + JI/2-1M - Ji 

\ rn=— 00 



(62) 

These two expressions for the energy density coincide for the flux $ = 1, which corresponds 
to a center vortex. Thus, the energy density as a function of the flux $ is continuous at the 
center vortex flux $ = 1. Furthermore, one also immediately recognizes that the relation 
ei($ = 0;x) = 62 ($ = 2; x) holds, which implies that the energy density of a Dirac string, 
$ = 2, is equal to the energy density of the perturbative vacuum, $ = 0. This is expected, 
since the Dirac string does not carry observable physical flux. Let us emphasize that this 
degeneracy in the energy density justifies the self-adjoint continuation chosen above. 
Since the energy eigenvalues depend only on the combination A(a)$ and summation is carried 
out over A(a) = ±1, the energy density is reflexion symmetric, i.e. stays invariant under the 
replacement $ — $. Furthermore, as can be read off from the explicit expressions 
the energy densities for the flux ranges < $ < 1 and 1 < $ < 2 are related by 



ei ($) = 62 (2 - $) . (63) 

We can therefore restrict ourselves to the flux interval < $ < 1. Subtracting from the 
energy density of the quark fluctuations in the presence of the vortex the corresponding 
energy density in the perturbative vacuum, one finds 
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Aei($;x) = ei(<l>;a:;) - ei(<l> = 0;a;) (64) 

oo oo ^ 



-oo 



oo oo 



^ y d/c y KdK^/K?T~k^ (.P^/^if^p) ~ Jl/^iKp)) (66) 



-oo 



The sums which occur here are identical to the sums which occur in the gluonic part of the 
energy density. In fact, using the results obtained for the gluonic part, see eq. and 
replacing $ by $/2, one finds for the first term 



1 $ / /<l>\^\ / $ 



Analogous calculations as in the gluon sector yield for the second term 



B= 2 



oo 

I dttV^^t^/^e- (/_,/,(t) - U/,{t)) . (68) 







The remaining integral has to be carried out numerically for given vortex flux $. Adding 
this contribution to the gluonic part, we obtain the total energy density of the fluctuating 
gauge field and of the fermions in the vortex background field. The result is shown in Fig. |21 



IV. DISCUSSION 



We have evaluated the one-loop effective energy density associated with an infinitely thin 
SU(2) Yang-Mills vortex, taking into account both gluonic as well as (massless) quark fluc- 
tuations, for arbitrary vortex flux $. As a function of the latter, center vortices (odd integer 
$) represent local minima of the effective energy - it is energetically favorable to keep mag- 
netic flux bundled in units corresponding to center elements of the gauge group as opposed 
to letting magnetic flux disassociate and propagate in lines of arbitrary flux through space. 
In this sense, the quantum fluctuations tend to stabilize center vortices with respect to local 
deformations. 
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FIG. 2: Total energy density of the fluctuating gauge field and of the fermions. 

Globally, the thin center vortices we consider are degenerate with the perturbative vacuum 
if one takes into account only gluonic fluctuations, cf. Fig. ^ This is natural, since these 
infinitely thin, odd integer $ configurations correspond to pure gauge transformations in 
pure Yang-Mills theory - they are just as invisible to gluonic fluctuations as Dirac strings 
(even integer $), since gluons are defined in the adjoint representation of the gauge group. 
This degeneracy is lifted when quark fluctuations are included. The energy density of a 
center vortex is raised above the energy density of the perturbative vacuum through the 
influence of the fermions (but it continues to represent a local minimum as a function of 
magnetic flux $). Only Dirac strings are then still degenerate with the perturbative vacuum. 
In a realistic vortex model of the QCD vacuum, energy densities of the type considered here 
are but one among a whole set of factors influencing vortex dynamics. For example, realistic 
vortices are not infinitely thin, but are smeared out to a finite transverse thickness, a case 
which has recently also been d.™.edQQ. AS Shown th.e,a«.tetMcU„e.ca. actually 
lead to a lowering of the effective energy density below the one of the perturbative vacuum. 
On the other hand, also vortex interactions have to be considered, cf. |22] for an investigation 
on the lattice. In the phenomenologically successful random vortex world-surface model 



IS In 



the energy associated with vortex curvature was found to play an important 
role - this contribution has not been investigated from first principles. In that model. 
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it is ultimately the entropy inherent in the random vortex world-surface ensemble which 
overcomes any energy penalties associated with the formation of center vortices, and thus 
leads to the formation of the percolating vortex vacuum and a comprehensive description of 
the principal nonperturbative phenomena characterizing the strong interaction. 
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APPENDIX A: CLOSURE OF C 

In this Appendix we discuss the closure of the operator C*, cf. (j2ip . We prove that functions 
/ G D{C) are absolute continuous close to the origin and fulfill lim^^QX'^ f{x) = 0, where 
S =\m- \{a)^\. To obtain D{C), we add to the domain of C, D{C) = C^{R>o), the limit 
points of Cauchy sequences in D{C), whose images under C are also Cauchy sequences. We 
consider such a Cauchy sequence fk and denote fki = fk — fi- We obtain 

oo oo 

J dx {xfli^ + -fki^) = J xdx fkiCfki < e , for A;, / > iV^ . 



Therefore, x~^l'^fk and x^^'^f'^ are Cauchy sequences in L^(]R>0) dx). Observe that also x^^'^fk 
and x^^'^Cfk are Cauchy sequences in L^(]R>o;dx). 

Next we prove that if x"'/^, and x^~^fk are Cauchy sequences {n G R) in L'^{Ii^Q,dx), 
then also x"/^"^/^/^ and x"/^"^/^/^ are Cauchy sequences. From that, it follows by induc- 
tion that x'^/^^"/^ and x^^^'^^^fk are Cauchy sequences for any < e ^ 1. We consider 
dx x^ fli^C fki < ^ for /c, Z > and obtain by partial integration 

oo 


We obtain another relation from Jq°° dx x'^~^ fki^/xC f^i < C , for k, I > N^: 

oo 

/ dx (x-'/'f'J + {6- ^^^^|^)x"-^/V..^) < i . 
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For 5 ^ Q, n e Q and vice versa, which can be achieved for any < 5 < 1 by adding an 
infinitesimal number to n, these two relations are linearly independent, as one can check by 
calculating the determinant. In these cases x^/^-^/^/^/ and x^^'^"^/'^ f'^i are Cauchy sequences 
in L^(]R>o, dx). 

In the next step, we prove that the functions in D{C) arc absolute continuous close to the 
origin. We consider an interval (0, a) and denote the limit of y^/jt by /, the limit of ^/xC fk 
by the hmit of a;~^/^+Vik by g and the limit of x^^^^'^^fk by g in L^(R>o, dx). We get 

a a 

J dx \h - x'/'-%' < a'-'^J dx \x-'/'^^h -g\'<^ 



that is 



fk ^ x'/^-^g = h 

in L^((0, a), dx). On the other hand, as gi e L^((0, a), dx) and therefore the integral 
Jq dx' x'^^'^^^g is well defined for all x e (0, a), we have 

2 

J dx' x"/'-^g -h < a'-^^ J dx' \g - x'^^+V^r < i 



and see that 



fk 



X 



x"/'-'g 



However, as — > /i we see that h is absolute continuous in (0, a) and -^t^ is absolute 
continuous in (a', a) for < a' < a. Taking the derivative of ^^27 we obtain 

i^J - 9j^,-{l-2e)gJ^^eL\{0,a),dx) 
as all functions on the right side are in L'^{{0, a), dx). Therefore 

a 

h 



dx 



X 



l-2t ' 



h(a) h(a') 
——k — hm — 



has to be finite and we arrive at 



lim h{a')a' 

a'— »0 



/V_/-l+2£ 



const < 00 
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for any < e ^ 1, i.e. 

x~^h{x) — > for X 

for all < 5 < 1. 
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